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Abstract. A simple two-species asymmetric exclusion model in one dimension with 
bulk and boundary exchanges of particles is investigated for the existence of spontaneous 
symmetry breaking. The model is a generalization of the 'bridge' model for which earlier 
studies have confirmed the existence of symmetry-broken phases, and the motivation here 
is to check the robustness of the observed symmetry breaking with respect to additional 
dynamical moves, in particular, the boundary exchange of the two species of particles. 
Our analysis, based on general considerations, mean-field approximation and numerical 
simulations, shows that the symmetry breaking in the bridge model is sustained for a range 
of values of the boundary exchange rate. Moreover, the mechanism through which symmetry 
is broken is similar to that in the bridge model. Our analysis allows us to plot the complete 
phase diagram of the model, demarcating regions of symmetric and symmetry- broken phases. 
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1. Introduction 



One-dimensional systems with short range interactions and a finite state space for the local 
variable, when driven far from equilibrium, often exhibit phenomena which are unexpected 
in equilibrium P, |2]. An example is that of spontaneous symmetry breaking which was 
observed in a driven system of two species of particles [H Hj. In this so-called bridge model, 
one considers a one- dimensional lattice with sites occupied by two species of hard core 
particles, referred to below as the positive and the negative particles. The positive particles 
move stochastically to the right, while the negative particles move stochastically to the left. 
At the left boundary site, positive particles may enter the lattice and negative particles may 
leave; at the right boundary site, negative particles may enter the lattice and positive particles 
may leave. The dynamics is symmetric with respect to "charge" conjugation combined with 
space inversion. 

At long times, the system reaches a nonequilibrium stationary state with non-zero 
particle currents. For low extraction rate of particles, the system is typically loaded with 
a majority species of particles, say, the positive particles, with larger current and higher 
bulk density than the minority species (the negative particles). As the system evolves, it 
fiips between this state and the one in which the current and the density inequalities are 
reversed. General considerations and Monte Carlo simulations showed that the average time 
between fiips grows exponentially with the system size, leading to spontaneous symmetry 
breaking (SSB) in the thermodynamic limit [3]. The occurrence of SSB for the bridge model 
was demonstrated within mean-field approximation [3l H] and also by rigorous results under 
specific conditions [SlEllTj. Over the years, such symmetry breaking has also been observed 
in many variants of the bridge model [H O [TOl [11], [121 [13 [H [IS] . It may be noted that, for a 
periodic system, the particle number is conserved and hence, SSB cannot occur. Therefore, 
SSB in the bridge model may be understood as a boundary-induced critical phenomenon. 

In this paper, we study the robustness of SSB in the bridge model with respect to 
additional dynamical moves at the boundary. It is known that far from equilibrium, 
boundary-induced critically phenomena generally depend sensitively on microscopic details 
of the boundary processes [16]. It is thus of interest to test the effect of more general 
boundary processes on the stationary state. To this end, we generalize the bridge model to 
allow for boundary exchange of particles of the two species. Our model is defined as follows. 

We consider a one- dimensional lattice of sites. Each site i of the lattice is occupied 
by either a positive (+) or a negative (— ) particle, or is left vacant, denoted by (a "hole"). 
The system evolves according to a stochastic Markovian dynamics. In an infinitesimal time 
dt, the following exchanges may take place at a pair of nearest-neighbor sites {i, i + in the 



bulk (1 < z < A^ - 1): 
(+).(0)m 



(0)j(+)j+i with probability dt, 
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with probability qdt. 

At the boundaries, particles enter or leave the lattice. During an infinitesimal time dt, the 
following events may take place at the left boundary site (z = 1): 

(0)1 (+)i with probability adt, 

(-)i (0)1 with probability /3dt, (2) 
(— )i (+)i with probability •ydt, 

while the following events may take place at the right boundary site (i = N): 

{0)n — ^ i~)N with probability adt, 

{+)n {0)n with probability /Sdt, (3) 
{+)n {~)n with probability •ydt. 

From the dynamical rules in Eqs. ([T]l3]), it is clear that the dynamics is symmetric with 
respect to charge conjugation (+ — ) combined with space inversion (left right). On 
setting the boundary exchange rate 7 to zero, we recover the bridge model studied in [3], H] 
which exhibits SSB, as discussed above. In these studies, detailed analysis of the case q = I 
was carried out. It was also argued that, despite some qualitative changes in the phase 
diagram as the bulk hopping rate q is varied, the phenomenon of SSB persists for q ^ 1. 
This is consistent with the notion that, in the bridge model, SSB is caused by boundary 
effects. We thus restrict our investigation of robustness of SSB in this paper to the case 
q = 1. Note that in our model, the dynamical moves are the most general ones consistent 
with the symmetry of the model and with the total asymmetry in the direction of motion of 
the particles. 

In our model, if the exchange rate 7 equals the particle injection rate a, a previous 
study has demonstrated that no symmetry breaking takes place [17] . Moreover, as explained 
below, if particle extraction occurs only by boundary exchange, i.e., if the particle extraction 
rate j3 equals zero, SSB disappears. Thus, the question arises as to how robust SSB is with 
respect to the boundary exchange process for general values of 7. 

In this work, we obtain the complete phase diagram of the model in the space of the three 
rates a,/3,7, and specify regions where SSB occurs. We base our analysis on a mean-field 
approximation and an exact analysis in certain parameter regimes. 

We find that for nonzero exchange rate 7 and for (3 0, the SSB in the bridge model 
is sustained so long as these rates are not too large. This statement is quantified later in 
the paper. Thus, SSB in the bridge model is indeed quite robust with respect to additional 
dynamical steps. Similar to the original model, we find that the average fiipping time between 
the two long-lived states in the symmetry-broken phases grows exponentially with the system 
size. Our results are corroborated by extensive Monte Carlo simulations of the model. 

The paper is organized as follows. In the following section, we discuss how our model 
is related to the single-species totally asymmetric simple exclusion process (TASEP), and 
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briefly summarize the known phase diagram for the TASEP with open boundary conditions, 
for use in later parts of the paper. In Section [31 we present exact results for the stationary 
density profiles and currents for two cases, namely, (i) when the particle extraction rate (3 
is zero, and (ii) when the boundary exchange rate 7 equals the particle injection rate a for 
non-zero extraction rate 13. In both these cases, the system exhibits symmetric phases only. 

In the absence of an exact solution for the stationary state measure for general 
values of the system parameters, in Section |H we make progress by applying a mean-field 
approximation to our model. In particular, we investigate the possibility of various symmetric 
and symmetry-broken phases in the stationary state. Combined with the results from the 
previous section, we obtain the complete phase diagram of the model, in which SSB occurs 
in a certain parameter regime. In Section [5], we report extensive Monte Carlo simulations 
for the particle density profiles in the stationary state, both for the symmetric and the 
symmetry-broken phases. We also provide numerical evidence for the exponential growth 
of the average flipping time with system size in the symmetry-broken phases. These results 
support the mean-field prediction of SSB in our model. In Section [6l we briefly describe, 
with the help of a toy model, the physical mechanism through which symmetry breaking 
occurs in our model for small values of the particle extraction rate. The paper ends with 
conclusions in Section [71 

2. Relation of our model to TASEP 

Our model generalizes the single-species totally asymmetric simple exclusion process 
(TASEP) to two species of particles. The TASEP is a paradigmatic model to study 
nonequilibrium driven systems [16]. On a periodic one-dimensional lattice, the model involves 
single species of particles, say, positive, moving stochastically round the lattice by exchanging 
with nearest-neighbor holes. Defining our model on a periodic lattice with just the bulk 
dynamics (Eq. ([T])) and no boundary dynamics, it is clear then that in our model, a positive 
particle, in its motion round the lattice, will not distinguish between a negative particle and 
a hole; also, a negative particle will not distinguish between a positive particle and a hole. 
Thus, for our model on a periodic lattice, the dynamics of the positive and the negative 
particles becomes that of two separate TASEP 's. Indeed, because of these rules, the two 
particle species (positive and negative) of the two TASEP 's behave microscopically entirely 
independently even though they share the same lattice. The difference between this single- 
lane model and two non-interacting TASEP's moving on two separate lattices appears only 
on a coarse-grained level of description through the constraint that the total particle density 
of positive and negative particles cannot exceed one. Note, however, that the dynamical 
moves of holes in our model on a periodic lattice do not become that of a TASEP particle. 
In fact, interpreting the positive particles as regular or first class particles and the negative 
particles as vacancies, the holes act as the so-called second class particles [H]. These second 
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class particles behave like the first class particles in exchanges with vacancies, while they act 
like the vacancies in exchanges with the first class particles. 

In our model on an open lattice with both bulk and boundary dynamics, the two 
TASEP's of the positive and the negative particles are not entirely independent since the 
two species interact microscopically at the boundary sites. It turns out that much of the 
behavior of the model can be deduced from the properties of the single-species TASEP with 
open boundaries. This process has been well studied in the past, and its phase diagram is 
exactly known [191 ISOl [21]. For use in the later parts of the paper, we briefly summarize 
below the phase digram of the TASEP in the limit of an infinite system. 

Let denote the injection rate of particles at the left end of the lattice, while stands 
for the extraction rate at the right end of the lattice (Here, the subscript s refers to single 
species.). The particle exchange rate in the bulk is set to 1. In the thermodynamic limit, 
the phase diagram comprises three phases. 

(i) A maximal current or power law phase for > 1/2, /3s > 1/2. In this phase, the particle 
density approaches from the boundaries to the bulk value of 1/2 as a power law, and the 
current is maximal {js = 1/4). 

(ii) A low-density phase for as < Ps, < 1/2- Here, the particle current is js = ^^(l — a^); 
the bulk density equals as{< 1/2), and is approached exponentially from the right boundary. 

(iii) A high-density phase for Ps < c^s^Ps < 1/2- Here, the particle current is js = /9s(l — Ps)', 
the bulk density equals 1—Ps{> 1/2), and is approached exponentially from the left boundary. 

The transition from the low-density and the high-density phase into the power law 
phase is continuous, while the transition from the high-density to the low-density phase is 
discontinuous, with phase coexistence on the line = /3s < 1/2 for which the density profile 
is linear. The linear density profile is a result of superposition of configurations in which 
a left-hand region of density as coexists with a right-hand region of density I — as with a 
shock between them [T9| [2T| [22] . 

3. The phase diagram: Exact results on specific planes 

In this section, we discuss two cases for which exact results for the stationary state current 
and density profile of particles can be obtained for our model. This is possible because, in 
these two cases, the dynamics of the system effectively becomes that of either the single- 
species TASEP or two decoupled TASEP's. The two cases are (i) the a — 7 plane (i.e., 
P = 0), and (ii) the plane = 7, with P non-zero. In both these cases, the system exhibits 
only symmetric phases in which the magnitude of the average positive current (denoted by 
j"*") equals that of the average negative current (denoted by j~). 
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3.1. The a — 7 plane with /9 = 



If the extraction rate (3 is 0, from the dynamical rules in Eqs. ([T]l3]), it is evident that holes 
cannot be injected into the system, and, thus, the system will have just the positive and 
the negative particles in the stationary state. If one now interprets the negative particles as 
"holes", the dynamics becomes that of a single-species TASEP with the positive particles 
being injected and extracted with the same rate 7; the system will thus exhibit the following 
symmetric phases, (i) For 7 > 1/2: Maximal current phase for the positive and the negative 
particles, with j"*" = j~ = 1/4. The density for both the particles equals 1/2 in the bulk, 
and is approached from the boundaries as a power law. In the spirit of [3], we refer to this 
phase as the power-law (pi) phase, (ii) For 7 < 1/2: Here, j"*" = j~ = 7(1 — 7), and the 
density profile is linear for both the species. This phase corresponds to the coexistence line 
= /5s < 1/2 of the single-species TASEP phase diagram. We refer to this phase as the 
coexistence phase. The schematic phase diagram for the a — 7 plane with /3 = is shown in 

Fig. m 
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Figure 1. The phase diagram on the a — 7 plane for /3 = 0, based on the exact analysis in 
Section [XTl 



3.2. The plane = 7 with f3 non-zero 

In this case, the dynamical rules in Eqs. ([UlSj) imply that both in the bulk and at the left 
(respectively, right) boundary, a positive (respectively, negative) particle does not distinguish 
between a negative (respectively, positive) particle and a hole. As a result, the two TASEP's 
of the positive and the negative particles are decoupled. The injection rate for both the 
particles is a, while the extraction rate for both is a + (3. The following symmetric phases 
are possible, (i) For a > 1/2: Power-law (pl) phase, with = j~ = 1/4. The bulk density 
equals 1/2 for both the particles and is approached from the two boundaries as a power law. 
(ii) For a < 1/2: Low-density in the bulk for both the positive and the negative particles. 
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with j"*" = j~ = a{l — a). The density of the positive (respectively, negative) particles 
decays exponentially to the bulk value a as one moves away from the right (respectively, 
left) boundary. Following we call this the low-density (Id) phase. The transition from 
the Id to the pi phase is continuous. The schematic phase diagram for the plane = 7 with 
non-zero /3 is shown in Fig. [2l 




Figure 2. The phase diagram on the a = 7 plane for /3 7^ 0, based on tlie exact analysis in 
Section [321 



4. The Full Phase Diagram: A Mean-field Study 

In this section, we study the full phase diagram of our model to identify the possible 
symmetric and symmetry-broken phases. In the absence of an exact solution for the 
stationary state for general values of the rates a,/?, 7, we pursue our study by employing a 
mean-field approximation. In this approximation, pair and higher-order correlation functions 
for particle occupations are approximated by products of average occupation numbers. 

Let us define two occupation numbers, and 6i, for each site i, where 
Tj (respectively, ft/) = 1 if site i has a positive (respectively, negative) particle and is zero 
otherwise. The hard-core constraint implies only one of Tj and 6i to be non-zero at a time. 
The occupation number of holes at site i is given by CTj = 1 — Tj — 6'j. In this notation, the 
magnitudes of the positive and the negative currents in the bulk {1 < i < N — 1), derived 
from Eq. ([1]) (with q = 1), are 

Ji|i+i = {r^(r,+l) + {tA+i) = {nil - T-,+i)), 

(4) 
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At the two boundaries, the currents (magnitudes only) are given from Eq. ([2]) and Eq. ([3]) 

by 

Jo+i = a{ai) + 7(^1) = «(1 - (ri)) - (a - 7)(^i), 

(5) 

Jn,n+i = (/5 + 7)(t-7v), 

jN+i,N = <^{(^n) + l{rN) = a(l - {On)) - (a - 7)(t"7v)- 

In the above, we have considered only the magnitudes of the current with the 
understanding that the positive current flows from left to right, while the negative current 
flows from right to left. Let pi (respectively, mj) denote the average density of the positive 
(respectively, negative) particles at site i. 

Pi = in), 

(6) 

rrii = {9i). 

These densities evolve in time according to 

dpi_ ^ ■+ _ .+ 
dt -^^'^+1' 

(7) 

drrii 

= 3i+i,i- 

We proceed by employing the mean-field approximation, in which the currents on the 
right hand side of the above equation are replaced by their mean-field values which are 
obtained from Eq. (j4j) and Eq. ([5]) as 

(8) 

3i+i,i =mi+i(l - mi), 
for the bulk (1 < i < iV - 1), and 

icti = a{l - pi) - {a - -i)mi, 
Jlfl ={P + l)mu 

(9) 

for the boundaries. 



8 



In the stationary state, the currents of the positive and the negative particles will be 
constant throughout the system, i.e., j^^^^i = and J^f i « ~ i~ ^- ^ result, for 

1 < z < — 1 , one has 

j+ =Pi{l -Pi+i), 

(10) 

j~ = mi+i(l - rrii), 
whereas, at the boundaries, one has 

j+ = a{l - pi) - (a - 7)mi = (/3 + -f)pN, 

(11) 

j~ = {p + 7)mi = a(l - mjv) - (a - 'j)pn- 

Following [3], we define for the positive species, the effective injection rate at the left 
boundary and the effective extraction rate l3~^ at the right boundary in the following way. 

/3+ = —. (13) 

Pn 



From Eq. (ITTI) . it follows that 



(14) 

/?+ = /? + 7. 

Similarly, one defines for the negative species, the effective injection rate a~ at the right 
boundary and the effective extraction rate (3~ at the left boundary in the following way. 

«- = -y (15) 

(1 - mN) 

= ^. (16) 
mi 

On using Eq. ffTTl) . we get 

J" 



(17) 

Thus, we have two TASEP's of the positive and the negative particles. The positive particles 
enter at the left boundary with rate hop through the bulk with rate 1, and exit at the 
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right boundary with rate jS^ . On the other hand, the negative particles enter at the right 
boundary with rate a~, hop through the bulk with rate 1, and exit at the left boundary with 
rate (3~ . The two TASEP's are coupled at the boundaries through the effective injection rates 
a^, which, for one species, depend on the density of the other. Using the phase diagram 
of the single-species TASEP, we discuss below the possible phases in our model, first the 
symmetric ones, followed by those with spontaneous symmetry breaking. In the following, 
we denote the bulk densities away from the boundaries by p and m for the positive and the 
negative particles, respectively. 



4--1- Symmetric Phases 

The two possible symmetric phases are the power-law (pi) phase and the low-density (Id) 
phase. One cannot have a high-density symmetric phase, since, in that case, the combined 
bulk density of particles would exceed one. In these symmetric phases, one has = 
(= j*, say) and therefore, from Eq. f|T^ and Eq. f|T7j) . one has 

, _ a(/9 + 7) ^ 

a = a = -— = a , say, (18) 

a + p 

while one already has l3~^ = j3~ in the definition of the model. 

• pi phase: 

Here, = 1/4, and far away from the boundaries, the densities are p = m = 1/2. The 
conditions for the occurrence of this phase are 

a' > l/2,/3 + 7 > 1/2. (19) 

Now, /3 + 7 is always greater than or equal to a^, so that there is a single condition 
defining this phase, namely, 

a' > ^ (pi phase). (20) 

• Id phase: 

In this phase, = a*(l — a*), and the bulk densities are given hj p = m = a*. This 
phase exists provided the following conditions are satisfied. 

a' < p + -f,a' <l/2. (21) 

The first condition is always satisfied, so that the single condition defining this phase is 

a' < 1/2 (Id phase). (22) 

From Eq. (1201) and Eq. (1221) . it follows that the surface in the {a, /3, 7) space that marks the 
transition between the two symmetric phases is given by 

„• = = 1. (23) 
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Solving for /3, one gets 

The phase transition from the Id to the pi phase is continuous, since, across the transition 
point, the current j'^ and its first derivative with respect to density are continuous, while the 
second derivative is not. 

4-2. Symmetry-broken Phases 

A symmetry-broken phase necessarily means unequal densities for the positive and the 
negative particles. Then, the different possible symmetry-broken phases are given by the 
six combinations of the single-species TASEP phases, summarized in Section [2], for the 
positive and the negative particles. Of these, the high-density/high-density phase and the 
high-density/power-law phase cannot occur, since for these combinations, the combined bulk 
densities of the positive and the negative particles would be greater than 1. It can be shown 
that the power-law/low-density phase does not exist within the mean-field approximation; 



the proof is outlined in [Appendix A[ Thus, the possible symmetry-broken phases are 
the high-density/low-density (hd/ld) phase and the low-density/low-density (Id/ld) phase. 
Below we investigate the possibility of occurrence of these phases within the mean-field 
approximation. We assume, without loss of generality, that the positive particles are in the 
majority in these phases. 

• hd/ld phase: Here, we have = (/? + 7)(1 — /3 — 7) and j~ = a~{l — a~). The bulk 
densities are given hj p = 1 — /5 — 7, m = a^. The conditions for the existence of this 
phase are 

+ 7, /3 + 7< 1/2, 

(25) 

a~ < /3 7, a" < 1/2. 

• Id/ld phase: 

Here, both the positive and the negative particles have unequal bulk densities of values 
smaller than 1/2. In this phase, we have = a~^{l — a~^) and j~ = a~{l — a~). The 
bulk densities are given hj p = ,m = a~ . With the positive particles in the majority, 
one has > a~ , and consequently, > j~ ■ The conditions for the existence of this 
phase are 

a+ < /3-F7,a+ < 1/2, 

(26) 

a~ < (3 + 'y,a~ < 1/2. 
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Comparing the conditions in Eqs. fl25H26l) . it follows that the transition surface between 
the hd/ld phase and the Id/ld phase is given by 

a+ = (3 + -f. (27) 

To plot the above surface, one has to express in terms of a, (3, and 7. In the hd/ld phase, 
one has, from Eq. (ITTI) . on substituting j+ = (/3 + 7)(1 — /3 — 7) and j~ = a~{l — a~), 

«" = ^Z^^^Zfl , (28) 

a-{l~a-) , / a-7\ (/3+7) (1-/3-7) 
a ~^ V/3+7 / a 

which gives a quadratic equation in a~ . Solving this equation, we get 

«- = - ^ V(l + - 4{«(/3 + 7) + 7(1 - /5 - 7)}, (29) 

where the root with the negative sign is taken to ensure that a~ < 1/2. The above equation 
gives the value of a~ throughout the hd/ld and the Id/ld phase. 

In order to get an expression for for the hd/ld phase, we substitute in Eq. 0141) the 
currents j"*" = (/5 + 7)(1 — /9 — 7) and j" = a~(l — a~) for this phase; we get 

(/3 + 7)(l-/3-7) 

= 7 r . 30) 

(/3+7) (1-/3-7) , f a-{l-a-) 



,/3+7, 

For the Id/ld phase, the expression for a+ follows from Eq. (1141) on substituting the currents 
for this phase, namely, j+ = a~^{l — a"*") and j~ = a~{l — a^); one gets 

+ a+(l -a+) 
= 7 r . (31) 

a+{l—a+) I / a— 7 \ a~{l—a~) 
a \/3+7 y Q 

The above equation leads to a quadratic equation in a"*", which is solved and the root with 
the negative sign, satisfying a"*" < 1/2, gives a+ in the Id/ld phase. 

From Eqs. (l30l) and ( 13T|) . it is clear that the expressions for for the hd/ld and the 
Id/ld phases become identical when the condition a"*" = /3 + 7 is satisfied; this condition thus 
gives the transition surface between the two phases, as already discussed above. In order 
to get an explicit expression for this surface in terms of a,/?, 7, one may proceed as follow. 
Inserting the expression for a~ from Eq. fl29l) into Eq. fl30l) . evaluating a+ and then equating 
to /5 + 7, we get a cubic equation in a in terms of (3 and 7. We solve this equation, and 
choose the particular root which gives as a function of a and 7 such that the condition, 
/3 + 7 < 1/2, is satisfied; we finally get the equation of the surface separating the hd/ld phase 
from the Id/ld phase. 

The existence of the Id/ld phase in our model is demonstrated within the mean-field 
approximation. In the original bridge model, whether this phase exists beyond the mean-field 
approximation has been a subject of some debate [23l [2^ [25]. 

It may be checked from Eqs. dTl), 1^, ([29]), ([30]) and that in the regions 
corresponding to the intersection of the hd/ld phase and the Id/ld phase with the (3 = 
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plane, one has = a" = = P~ = 7. Thus, there is no physical significance of the 
intersection of the transition surface, a"'" = /3 + 7, with the f3 = plane. This is because, 
on either side of the intersection curve = 7, the injection rate and the extraction rate of 
both the positive and the negative particles equal 7. This is consistent with the observations 
in Section 13.11 and consequently, one has only symmetric phases on the /3 = plane. 

4-3. Transition between the Id and the Id/ld phases 

From Eq. ( fT4l) and Eq. ( fT71) . on substituting the current for the Id/ld phase, namely, 
j"*" = a"^(l — a"*") and j~ = a~(l — a~), we get 

a"*" = 1 a'^il — a~^) — —a~(l — a~), 

a ^ ' (3* ^ ' 



(32) 



_ 1 1 _ _ 
a =1 — —a^(\ — a^) a (1 — a ), 



where 



r ^ (33) 

a — 7 

To solve Eq. flH21) . we follow the procedure in [1], and define 

^ = + = (34) 

On taking the difference of the two equations in (132|) . we get 

D={'^^i\D{\-S). (35) 



D 



Now, since we are dealing with a symmetry-broken phase, D 7^ 0. Thus, one has 

S=l-^. (36) 
Summing the two equations in ( !32!) . we get 

At the transition to the Id phase, the densities become equal so that D = 0. Thus, either 
5 = 2, which is excluded, since this means that the sum of the bulk densities of particles is 
greater than 1, or, that 

(38) 

The above equation, combined with Eq. fl5B]) . gives the following expression for the rates 
a"'", a~ on the transition surface between the Id/ld phase and the Id phase. 

a- = a- = ^ = ifl-^V (39) 



a + /9* 2 V a - /5* 
13 



Solving the above equation for /3, we get 



2 1 - a 



In the hmit a — > 1, we get 



/S^i^. (41) 

The intersection of the transition surface in Eq. (HOj) with the /3 = plane is along the curve 
2(27^ - 7) 

Note, however, that there is no physical significance of the above intersection curve. 
This is because one can check that on either side of this intersection curve, the injection rate 
and the extraction rate of both the positive and the negative particles are equal to 7. This 
observation is consistent with that in Section [3TT1 and consequently, one has only symmetric 
phases on the /3 = plane. 

4.4- Summary of the phase diagram 

We give below the equations of the three surfaces separating the various phases within the 
mean-field theory. 

• Surface that separates the Id phase from the pi phase: 

^="'^-^^>. (43) 

• Surface that separates the Id phase from the Id/ld phase: 

1 -47 
B = for a = 1. 

' 3 

(44) 

-3a2 + 5a7 + (VOa^ - 4a + 4) (a - 7) - 27 
~ 2(1 - a) 

for a 7^ 1. 

• Surface that separates the Id/ld phase from the hd/ld phase: 

a+ = /3 + 7, (45) 
where a"*" may be found from Eqs. ( l29l) and ( l30ll . 

Note that the two transition surfaces, given by Eqs. ( l44l) and ( H5l) . coincide on the /3 = 
plane. This can be seen in the schematic phase diagram in the a — /5 — 7 space, given in 
Fig. [3l In this phase diagram, we show only the regions occupied by the symmetry-broken 
phases, as predicted by the mean-field approximation to our model. In Fig. HI we show the 
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phase diagram on the — 7 plane for (3 = 0.01, obtained from our mean-field analysis. In 
Fig. O we show the mean-field phase diagram on the a — (3 plane for six values of 7. 

The transition from the hd/ld phase to the Id/ld phase is discontinuous, while that 
from the Id/ld phase to the Id phase is continuous. This is because, in the former case, 
the first derivative of the average current with respect to density is discontinuous across 
the transition, while, in the latter case, the second derivative of the current with respect to 
density changes discontinuously across the transition. 




Figure 3. Schematic pliasc diagram in the a — /3 — 7 space, showing only the regions 
occupied by the symmetry-broken phases, as predicted by the mean-field approximation. 
Within the region enclosed by the inner surface, one has the hd/ld phase, while, within the 
region between the inner surface and the outer shaded one, one has the Id/ld phase. The 
two surfaces intersect on the /3 = plane. As explained in the text, the hd/ld phase extends 
up to but not including the /3 = plane. 



5. Monte Carlo simulations 

To check the general features of the phase diagram predicted by the mean-field theory, 
and in particular, the existence of the symmetry-broken phases, we carried out extensive 
Monte Carlo (MC) simulations of our model. For a lattice of size and given values 
of the parameters a, (3 and 7, we followed a random sequential update scheme for the 
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Figure 4. The phase diagram on the a — 7 plane for /? = 0.01, based on the mean-field 
analysis in Section 2) The Id/ld phase occupies a very narrow region, and hence, appears 
as a line on the scale of the figure. 



configurations, with A^-|- 1 updates constituting one Monte Carlo Step (MCS). We typically 
ran the simulation for sufficiently long time (~ A^^ MCS) to ensure that the system relaxes 
to stationarity in this time, after which we started making measurements for the average 
density profiles and the currents for the two species of particles. 

In the symmetry-broken phases, in order to get the average density profiles, the 
occupation of each site was averaged over several runs of the simulation in the stationary 
state. However, one has to be careful that the running time of the simulations does not 
exceed the flipping time t{N) between the two states of the symmetry-broken phases. An 
estimate of t{N) may be made following a procedure explained later in this section. In the 
symmetric phases, no such restriction on the running time is necessary so that one runs the 
simulation long enough to reduce fluctuations in the measured density profiles. Following 
[1], the currents are measured as 

.± = ^ (46) 

where N± is the total number of positive (negative) particles which have moved in N^t 
MCS/site. 

The Monte Carlo density profiles for the symmetry-broken phases are shown in Figs. [6] 
and [7] corresponding to the hd/ld phase and the Id/ld phase, respectively, while those for the 
symmetric phases are shown in Figs. [Eland [9] corresponding to the Id phase and the pi phase, 
respectively. From the figures, it can be seen that the density profiles are fiat in the bulk, 
with some structures near the boundaries. In the hd/ld phase, the density of the positive 
particles is higher than 1/2, while the density of the negative particles is lower than 1/2; in 
the Id/ld phase, both densities are smaller than 1/2 and unequal in magnitude (shown in 
the blowup of the bulk density profiles). In the symmetric phases, on the other hand, the 
densities of both particle types are equal, taking a value which is either smaller than 1/2 (Id 
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Figure 5. The figure shows the mean- field phase diagram on the a — (3 plane for six values 
of 7. The phase boundaries are given by Eqs. (|43ll45p . The Id/ld phase occupies a very 
narrow region which appears as a line on the scale of the plots. 



phase), or, equal to 1/2 (pl phase). 

In order to illustrate symmetry breaking in our model, we present in Fig. [10] the time 
evolution of the current difference —j~ for a typical run of the MC simulation in the hd/ld 
phase. From the figure, it is evident that, excepting for short time intervals during which 
flips take place, the system is loaded predominantly with either the positive or the negative 
particles, alternating between the two as time progresses. We now proceed to show that 
the average time t{N) between successive flips of the current difference grows exponentially 
with the system size A^. In order to estimate t{N) numerically, we averaged the current 
difference over many runs, starting from the configuration where all sites are occupied by 
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1 256 512 

Site i 



Figure 6. Monte Carlo density profiles for the hd/ld phase. Legends: p ^ positive particles, 
m negative particles, h — > holes. Here, a = 1.0, f3 = 0.05, 7 = 0.1, N = 512. 




I ' 1 

1 256 512 




100 275 450 

Site i 

Figure 7. (a) Monte Carlo density profiles for the Id/ld phase. Legends: p positive 
particles, m negative particles, h — > holes. Here, a = 1.0, /? = 0.1984,7 = 0.1, = 512. 
(b) A blowup of the particle density profiles in the bulk, showing unequal values of the 
average density for the positive and the negative particles. 

one species of particles, say, the positive particles only [3]. This averaged quantity decays in 
time because of flips in current difference as a function of time. At large time t, in a system 
of size N, the average current difference decays as exp[— t/r(A^)], which gives the flipping 
time t{N). 

The results for t{N), extracted from MC simulations for various and given values of 
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Figures, (a) Monte Carlo density profiles for the Id phase. Legends: p ^ positive particles; 
m negative particles, h holes. Here, a = 1.0, (3 = 0.3, 7 = 0.1, N = 512. (b) A blowup 
of the particle density profiles in the bulk. Here, we plot pi and m^-i to demonstrate that 
charge conjugation combined with space inversion symmetry is not broken in this phase. 




Figure 9. Monte Carlo density profiles for the pi phase. Legends: p positive particles, 
m negative particles, h holes. Here, a — 1.0, f3 = 0.9, 7 — 0.1, N ~ 512. 



a,/3,7, is shown in Fig. [TTl It is readily seen that the plot asymptotically becomes hnear, 
indicating that the time scale t{N) grows as an exponential in the system size N for large 
N. This fact, combined with our observation for fiat density profiles in the bulk for the 
symmetry-broken phases, leads us to conclude that in the relevant parameter regimes (see 
Fig. [3]), our model exhibits spontaneous symmetry breaking. 
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Figure 10. Time evolution of the current difference in a typical run of the Monte Carlo 
simulation in the hd/ld phase. Here, a = 1.0,(3 = 0.05,7 = 0.1, = 256. Each point 
represents an average of the current difference over 1000 Monte Carlo sweeps. 
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Figure 11. Average flipping time t{N) as a function of system size N. Here, a — 1.0, /? = 
0.05,7 = 0.1. The points are obtained from Monte Carlo simulations. The dashed line is to 
indicate that the plot asymptotically becomes linear. 

In passing, we note that for fixed and finite A^, in the two symmetry-broken phases, 
as time passes, the average density profiles and currents for the particles would appear 
progressively symmetric due to repeated flips between the two symmetry-related states. 
Hence, while identifying symmetry-broken phases, it would be more appropriate to look at 
symmetric combinations of currents or densities, e.g., the sum and the absolute difference 
of currents Figure [12] shows + J~)/2 and |j+ — J~|/2 as a function of /3 for 

a = 1.0,7 = O-l- The points are obtained from MC simulations of a system of size 320, 
while the continuous lines are the mean-field predictions of Section HI One finds a fairly good 
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agreement for the sum, although not for the difference. 
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Figure 12. (a) Currents + j~)/2 and — j~\/2 as a function of P for a = 1.0,7 = 
0.1, = 320. The continuous lines represent the mean-field results, while points correspond 
to Monte Carlo simulation results, (b) Blowup of the dotted box in (a), showing the mean- 
field behavior of (j+ + j^)/2 in the transition region from the hd/ld to the Id/ld phase. 



6. Mechanism for symmetry breaking 

In this section, we briefly discuss the utility of the toy model of [H [1] in explaining the 
occurrence of spontaneous symmetry breaking in our model. Specifically, the model helps 
us to understand how, starting from the symmetric phase on the f3 = plane for 7 small 
(< 1/2), turning on arbitrarily small /? results in symmetry-broken phases. The toy model 
was initially devised to gain insight into the flipping process between the two states in the 
symmetry-broken phases of the original bridge model which has 7 = 0. The crucial step 
was to identify that, for small particle extraction rate (3, typical configurations on a lattice 
of size N are those composed of three blocks: a left block with an integer number j of the 
negative particles, a right block with an integer number k of the positive particles, and a 
central block with N — j — k holes. Other configurations in which holes are present inside 
the particle segments do not play a role in determining the long-time behavior of the system 
and may be neglected. On the time scale set by the rate (3, the dynamics of the system will 
involve transitions between various three-block configurations. 
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In our model, the extraction rate equals /3 + 7. Now, consider the case of small /3 + 7, 
with /3 <C 7. Typical configurations will again have three blocks, albeit with important 
differences from those in the original bridge model. On the time scale of most of the 
positive particles entering the lattice through exchanging the negative particles at the left 
end will have enough time to move away from the left end, so that the left block will have 
predominantly negative particles with the positive particles clustering towards the right end 
of the block. The crucial point to remember is that the left block will basically have no 
holes. Similarly, the right block will have predominantly positive particles towards the right 
end with some negative particles progressively clustered towards the left end of the right 
block. Again, the right block has basically no holes. All the holes are trapped between the 
left and the right blocks. Thus, the central block will have mostly holes with a few particles 
of both species here and there between the holes. A typical configuration will thus look like 

+ + + + -- + + - + + + +0 + 00- 0000 + - + 000- + + -- + 

H 1 1- + + + + +. The block of the negative particles at the extreme left end and the 

block of the positive particles at the extreme right end are long-lived in the limit of small 
/3 + 7. In the limit of small (3, the dynamics of the system will involve transitions between 
such three-block configurations. 

Now that we have identified the left block of length j with predominantly negative 
particles and the right block of length k with predominantly positive particles, the arguments 
given in the context of the toy model equally apply to the problem at hand. Thus, the 
dynamics restricted to {j, k) configurations will be that of a random walker in the first 
quadrant of the {j,k) space within the triangle with corners at (0, 0), (A^, 0), and (0,A^). 
Starting with a system with only positive particles (i.e., from the point (0,A^)), typical 
trajectories of the random walker will be biased towards the k axis [5]; see also [1]. In order 
to fiip to a system with only the negative particles, the trajectory will have to perform an 
atypical walk against the bias, starting from the point (0, A^) and ending on the j axis without 
touching the k axis. The probability of such a walk was calculated to be exponentially small 
in the system size [5]. As a result, the fiipping time t{N) diverges exponentially with the 
system size, as confirmed by the results of our simulation (Fig. [TT]) . Thus, the toy model 
of a biased random walker serves as a guide to understand intuitively the occurrence of 
spontaneous symmetry breaking in our model in the limit of small (3 and 7, with /? -C 7. 

7. Conclusions 

In this paper, we have revisited the problem of spontaneous symmetry breaking under 
nonequilibrium conditions in a one-dimensional model system with local dynamics and a 
finite state space for the local variables of the system. In this so-called bridge model, where 
two species of hard core particles are driven in opposite directions on an open lattice, earlier 
studies have shown the occurrence of SSB in the limit of small extraction rate f3 of particles 
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[3], H] . In this study, we have allowed for boundary exchange of particles of one species into 
another with rate 7. This is the most general dynamical move at the boundaries which is 
consistent with the symmetry of the model and with the total asymmetry in the direction 
of motion of the particles. In this modified model, if the rate /3 is zero, exact results predict 
only symmetric phases. On the other hand, if the rate (5 is large, one expects, on physical 
grounds, only symmetric phases. Also, if the exchange rate 7 equals the particle injection 
rate a, there is no SSB, as confirmed in [T7j. Thus, one is left to wonder about the occurrence 
of SSB for small values of (3 and non-zero 7. 

Here, we reported exact as well as mean-field results for the complete phase diagram of 
the modified bridge model in the a — /3 — 7 space, showing regions for the symmetric and the 
symmetry-broken phases. Our results confirm the existence of SSB for non-zero (3 and 7, 
provided both the rates are not too large, as have been quantified in the paper. Similar to the 
original bridge model, in the symmetry-broken phases, the system resides in one of two long- 
lived states, whose average lifetime grows exponentially with the system size. Our results are 
supported by extensive Monte Carlo simulations of the model. All these observations lead 
us to conclude that the symmetry breaking in the original bridge model is quite insensitive 
to additional dynamical moves, like those allowing for particle exchange at the boundaries, 
for a range of values of the exchange rate. It is left as an open problem to obtain the exact 
solution of the stationary state for the entire range of values of the parameters defining the 
model, for example, by including the boundary exchange rate 7 in the studies pursued in 
[6],[7j, and also, to study the model in higher dimensions. 
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Appendix A. Proof of non-existence of the power-law/low-density (pl/ld) phase 
within the mean-field approximation 

The conditions for the existence of the pl/ld phase (with positive particles in the pi phase) 
are 

a+ > 1/2, /3 + 7>l/2, a-</5 + 7, a" < 1/2. (A.l) 
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Correspondingly, the particle currents satisfy = 1/4 and j < 1/4. From Eq. (fTT]) . it 
then follows that 

j+ = a{l - pi) - (a - 7)mi = {(3 + -f)pN = 1/4, (A.2) 

and 

J- = a{l - mjv) - (a - 7W = (/5 + 7)"^i < 1/4. (A.3) 

In order to prove the non-existence of the pl/ld phase, we will use the above expressions 
for the currents as well as the facts that the negative particles are in the low-density phase 
so that < mj\f < 1/2 and that they have a flat density profile near the right boundary so 
that j" = mN{l — rriN). We get an expression for ttin which in turn yields an expression for 
mi. Now, since the positive particles are in the power-law phase, pi > 1/2, which, on using 
Eq. flA.2l) . puts a bound on mi. We will show that the derived expression for mi fails to 
satisfy this bound in the relevant parameter regime. 

From Eq. (lA.2p . we get 

Since the negative particles are in the low-density phase for which the density profile is fiat at 
the entrance end, one has j~ = m7v(l — mp^). Then, from Eq. (lA.Sp . on using the expression 
for ptv in Eq. (1A.4I) . we get 

(tt — 7) 

mAr(l - m^v) = a(l - mAr) - — — — -, (A. 5) 

4(p + 7) 

which gives a quadratic equation in rriN, and therefore. 



Here we have taken the negative root for ttin- This is because the negative particles are in 
the low-density phase, and hence, we should have < m^v < 1/2. The negative root satisfies 
these bounds, provided 

and 



The inequality in Eq. ( 1A.9I) may be rearranged to give 

/3 < (A.10) 
2a — i 
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so that, in order to have a finite /5 > 0, we need to have either a > 1/2, < 7 < 1/2, or 
< a < 1/2, 7 > 1/2. Also, we have, from Eq. f[0]) . 

/3 + 7>^- (A.ll) 
Let us introduce a new variable x by the following equation, 

(/9 + 7) 

For X > (when a > 1/2 and < 7 < 1/2), Eqs. flAlTl) . flAlSl) . fICT and flATTl) give the 
following bounds on x, 

2a-l<x<2a. (A.13) 

On the other hand, for x < (when < a < 1/2 and 7 > 1/2), Eqs. flXT]) . f[OD . (lAlQl) 
and (lA.lip allow for x to lie in the following range: 

2a-l<x<0. (A.14) 

Next, we derive an expression for nii from Eq. flA.SP by substituting j~ = miy{l — 171^), 
with niN given by Eq. (1A.6I) . We get 



1 



mi 



4(/? + 7) 



1 - a 



(l-a)2 + 



(q-7) 
(/5 + 7) 



On the other hand, Eq. (1A.2P yields 



1 



mi 



(a — 7) ^ ^ 4(^ — 7) 

Since the positive particles are in the power-law phase, we have pi > 1/2, i.e., I — pi < 1/2. 
Combined with the last equation, this gives 



(A.15) 



(A.16) 



mi < 



2a 



and 



mi > 



4(a - 7) 



2a 



for q; > 1/2 and 7 < 1/2, 



for a < 1/2 and 7 > 1/2, 



(A. 17) 



(A.18) 



4(7 - a) 

with mi given in Eq. ( 1A.15I) . In terms of the variable x, it follows from the above inequalities 
that, for positive x with 2a; — 1 < x < 2a, we must have 

2" 



X 



a 



a] 



X 



(2a- 1) < 0, 



(A.19) 



while, for negative x, the above inequality has to be satisfied for 2a — 1 < x < 0. It can 
be checked that in either case, the inequality has no solution, implying that the pl/ld phase 
does not exist in our model within the mean-field theory. 
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